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Выполнен расчет эволюции пластической деформации у вершины трещины в кристалле а-Ее в условиях 
плоской деформации (моды Г и П). Получены временные распределения пластической деформации, эффЕк- 
тивного сдвигового напряжения, коэффициента интенсивности напряжения. Проведено сравнение результа- 
тов расчетов для трещин с затупленной и острой вершинами. 
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Введение. Оценивание влияния формы вершины трещины на характеристики разрушения явля- 
ется актуальной задачей физики прочности и механики разрушения. Выбор модели формы хруп- 
кой трещины существенно влияет на результаты расчета распределения упругого напряжения и 
деформации в окрестности ее вершины. В простейшем случае трещины-разреза, когда пренебре- 
гают радиусом кривизны вершины р-> 0, эти распределения можно представить в виде асимпто- 


тического ряда по степеням расстояния „““'/? (к>0) от вершины до заданной точки [1, 2]. Пре- 
небрежение членами ряда с К <0 обусловлено условием ограниченности величин перемещений 
и энергии деформирования у вершины трещины. Расчеты показывают, что распределения напря- 
жений и деформаций содержат лишь особенности в вершине 0(+"°) при А =0, которые могут 
исчезнуть при некотором распределении внешней нагрузки [3]. Кроме детального изучения глав- 
ного члена асимптотики (А =0), содержащего особенность, в последние годы внимание исследо- 


вателей привлекает учет несингулярных слагаемых упругого поля вблизи вершины трещины- 
разреза (К >1). В работах [4, 5] подробно исследован второй член асимптотики упругого напря- 
жения (Т-напряжение), который в случае трещины-разреза не зависит от г. В этом случае 
Т-напряжение входит лишь в компоненту тензора напряжения, соответствующую растяжению или 
сжатию вдоль линии трещины. Расчеты показали, что учет несингулярной составляющей заметно 
изменяет условия роста трещины [4, 5]. 

Впервые оценка распределения упругого поля в полярных координатах ги ф у вершины 
хрупкой трещины с учетом ненулевой величины радиуса вершины р прямолинейной трещины 
получена в исследовании [6]. Однако позже была получена более точная формула для распреде- 
ления компонент тензора напряжений у вершины [7]: 
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В формуле (1), аналогично исследованию [7], К, и Ки - коэффициент интенсивности 
напряжения (КИН) эквивалентной трещины, т.е. острой трещины (р -> 0) стой же длиной и ори- 


ентацией, как затупленная трещина. Сравнение результатов [6, 7] с результатами расчета мето- 
дом конечных элементов показало, что с увеличением доли моды П во внешней нагрузке точ- 
ность аналитических формул для компонент тензора напряжений [7] выше. 

Для хрупкой трещины эллиптической формы при смешанном нагружении условие катаст- 
рофического роста трещины зависит не от формы трещины, а лишь от радиуса ее вершины р [8], 


который можно оценить как 
р=16,217,/о‹. (2) 

Для а-Ре удельная поверхностная энергия 7,=1,42 Дж/м? и о.=31,7 ГПа и оценка (2) дает 
р=0,726 нм (около 3 постоянных решетки). Отсюда следует, что хрупкую трещину в приближении 
механики сплошной среды, нужно рассматривать как разрез с р->0 [8]. 

В последние десятилетия обозначился значительный прогресс в создании дислокационных 
моделей для пластической деформации у вершины трещины [9, 10]. Расчеты взаимодействия 
дислокации с трещиной-разрезом обнаружили явление экранирования дислокацией ее упругого 
поля, которое характеризуется снижением КИН по сравнению с КИН трещины без дислокаций. 
Позже эти расчеты были обобщены для случая клинообразного выреза [11, 12], хотя в этих рабо- 
тах получена корневая особенность распределения напряжения у вершины трещины, соответст- 
вующая особенности трещины-разреза. В работе [13] исследовано взаимодействие краевой дис- 
локации с трещиной в форме сплюснутого эллипса. С использованием перехода к трещине- 
разрезу согласно исследованию [7] получены зависимости КИН от взаимного расположения тре- 
щины и дислокации и обнаружены положения дислокации, соответствующие максимальным эф- 
фектам экранирования (мода Г) и антиэкранирования (мода 11). В работе [14] исследованы зави- 
симости КИН модели трещины, состоящей из двух лунок, от положения винтовой и краевой дис- 
локаций относительно вершины. Здесь также получена корневая особенность в распределении 
напряжения в вершине трещины в форме двух лунок. Следует упомянуть также об оценках вели- 
чины взаимодействия дислокации с трещиной, в вершине которой находится зона сцепления 
[15, 16]. В этих работах показано, что если модель трещины с особенностью по напряжениям 
(высокая величина сцепления) в вершине предсказывает экранирование или антиэкранирование 
в зависимости от знака вектора Бюргерса дислокации, то при низкой величине сцепления - все- 
гда экранирование вне зависимости от знака вектора Бюргерса дислокации. 

Расчет равновесной формы дислокационной пластической зоны у вершины трещины- 
разреза (см., например, [17]) показал, что компонента напряжения оо на линии продолжения 


трещины достигает максимума на некотором расстоянии от ее вершины. А дислокации, испущен- 
ные из вершины, более эффективно экранируют упругое напряжение трещины, чем дислокацион- 
ные диполи, расположенные вблизи вершины [18]. Другой особенностью распределения дислока- 
ций в пластической зоне является подтвержденное расчетами существование свободной от дис- 
локаций области, прилегающей непосредственно к вершине трещины. 

Пластическая деформация у вершины трещины в нагруженном образце существенно зату- 
пляет трещину, а это влияет на условия зарождения дислокационных петель из нее [19]. Форма 
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вершины, в свою очередь, определяет распределение экранирующих дислокаций в пластической 
зоне. Взаимное влияние формы вершины трещины и распределение пластической деформации в 
ее окрестности существенно меняет оценку (2). В частности, взаимодействие трещины и пласти- 
ческой деформации у ее вершины определяет условие вязко-хрупкого перехода. 

Цель наших исследований —сравнение эволюции пластической деформации у вершины 
трещины в ОЦК кристалле при учете затупленности вершины трещины, условий нагружения и 
систем дислокационного скольжения при пластической деформации. В дальнейшем ограничимся 
расчетами для кристалла а-Ре (моды Ги П]) в условиях плоской деформации. 

Постановка задачи и метод решения. Рассмотрим плоские задачи об эволюции пластической 
деформации у вершины трещины длиной 2/, расположенной в плоскостях скола {100} и {110} 
бесконечного ОЦК кристалла с постоянной решетки а (рис.1). 


Рис.1. Кристаллографическая схема плоскостей 
легкого скольжения у вершины трещины 


В первом случае выберем направление линии фронта трещины <001> и систему легкого 
скольжения <111> {110}. В ней пластическая деформация осуществляется путем термоактивиро- 
ванного движения смешанных дислокаций с вектором Бюргерса Б=а/2<111>. Для трещины в 


плоскости {110} направления линий фронтов <111> и <110>, а системы легкого скольжения 
<111> {110} и <111> {112} соответственно (линии дислокаций параллельны фронтам трещины). 
Для этих систем скольжения пластическая деформация происходит с помощью тех же дислока- 
ций, как и выше, но во втором случае дислокации чисто краевые. В расчетах учитываются только 
краевые компоненты смешанных дислокаций. 

К плоскостям кристалла у=-+® приложены однородные напряжения растяжения 


с„(=о.(!) (мода Г) и сдвига ‹„(1)=т.(!) (мода П), монотонно возрастающие до некоторых 
значений ок, тик, достаточных для пластического деформирования кристалла, но недостаточ- 
ных для роста трещины. После достижения нагрузками о.(!) и т.(!) своих максимальных значе- 


ний расчет проводился в режиме релаксации напряжения у вершины трещины, которому соответ- 
ствуют постоянные внешние напряжения. В данном расчете принята пропорциональность двух 
типов внешних нагрузок с.(!) = Ат.(!) ‚ постоянная К является параметром задачи. 

В кристалле равномерно распределены источники дислокаций, испускающие прямоуголь- 


ные петли, лежащие в плоскостях легкого скольжения. В расчете учитываются только отрезки 
петель, перпендикулярные плоскости Оху ‚ И их краевые компоненты. Последнее упрощение свя- 
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зано с видом приложенной нагрузки, действующей только на данные компоненты дислокаций. За 
счет концентрации напряжения у вершины трещины, релаксирующего путем пластического де- 
формирования, в этой области достигаются значительные плотности эффективных дислокаций 
(избыточные дислокации одного знака среди дислокаций с параллельными векторами Бюргерса) 
и, следовательно, возникает упругое поле, существенно влияющее на эволюцию пластической 
деформации. 

Предположим, что пластические зоны, образующиеся в кристалле, имеют в плоскости 
Оху линейные размеры настолько малые по сравнению с /, что правомерно описание пластиче- 


ской деформации у вершины полубесконечной трещины. Плоскости {110} и {112}, пересекаясь с 
плоскостью 0ху, образуют на ней два семейства линий скольжения, направления которых заданы 


векторами &, (^=1, 2). 


Скорость пластического деформирования, обусловленного движением дислокаций у вер- 
шины трещины, дается формулой [20]: 


ме хр _ оо лы 
Я КвТо(г,Е) 


где И, - энергия активации; д, — постоянная Больцмана; Т. - температура; ви т, — постоян- 


яепо«(,!), (3) 


ные; о“ — эффективное напряжение сдвига в А-й плоскости легкого скольжения [20] 





01,0) = о\(",/)— а,(", Ото" (т, при |в“(». 1) > [6.(", 0) 


> 


с.(",Ю) =0 при |с*(",0)| <|6.(",0]. (9) 

В формуле (4) сдвиговое напряжение в плоскостях легкого скольжения 
(1,1) = о (Г, + (Р.И, (5) 
в,(,!) =со+ с (г,!) = (6) 


напряжение, препятствующее пластическому сдвигу за счет трения решетки сх и локального уп- 
рочнения материала с’, ‚ вычисляемого по формуле: 


б/ р] 


т 
, (7) 
Где о. И т - постоянные. 

В формуле (5) с+(и,Е) — сдвиговая компонента тензоров напряжения (1) для А-й плоскости 
легкого скольжения, которую в ряде расчетов дополняли вторым членом асимптотики 
(о, =о,+ Тсозф ‚ Оф=О+ ТзтФ ‚ оъ=т-ьы Тзшфсозф ), а также силами сцепления берегов 
трещины у ее вершины; с„(х,Г) — дальнодействующее упругое напряжение, создаваемое дисло- 
кациями одного знака в пластической зоне 

2 
о#(г,1) = х Гв(2)Ари2 ва,  г=хЧЬ, (8) 
ы В» 


где Р, - часть пластической зоны, образовавшейся в верхней полуплоскости в результате 
скольжения дислокаций вдоль &,. 
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В формуле (8) Ар, (2',Г) — плотность эффективных дислокаций, связанных с деформацией 
= (Р,Г) соотношением: 
Ари = ЕНи,о, (9) 
Ь ЧЕ, 

а 5,(2') определяет напряжение, создаваемое дислокацией в упругой плоскости с полубесконеч- 

ным разрезом, и задается выражением: 
би(т,"") =’, (г,"")+ о", (,"'), (10) 
где в'к(и,"') — компоненты тензора, которые характеризуют собственное поле пары дислокаций, 


а компоненты тензора "к (и,"') описывают поле напряжений полубесконечной трещины- 
разреза, нагруженной на ее берегах усилиями, равными по величине и противоположными 
по знаку усилиями, создаваемыми в сплошном кристалле напряжениями в'к(и,"') на месте 
трещины. 

Расчеты показали, что учет формы трещины не вносит существенных изменений в о%(г,г) 
по сравнению с трещиной-разрезом на расстояниях, больших радиуса вершины (’>р) или раз- 
мера зоны сцепления. 

Уравнения (9)-(17) образуют систему, из которой при начальных 

вк(гЕ=0)=0, о.(гё=0)=0 (11) 
и граничных условиях 
б(х<0,у=0,0 =0 (12) 
определяются =“(г,г), о“(г,г). 

Данную систему уравнений решали численно с изменяющимся в зависимости от вычисли- 
тельной ситуации шагом интегрирования Д;,‚. Автоматический выбор шага по времени учитывал 
ограничение на максимальную скорость деформации в формуле (3) плах (&(7,1)) < 0,1 ст. При 
этом ограничении остается справедливой формула (3), которая предполагает термоактивирован- 
ный механизм движения дислокаций. Метод решения системы (3)-(12) аналогичен методу в рабо- 
те [20]. В расчетах задавали не с» И т», а соответствующее значение максимального КИН 
Кг” и Ки. 

При расчете эволюции пластической деформации в кристалле о_Ее выбраны следующие 
значения постоянных: 21/=10° м, К" =0,2 МПа-м"?, Т.=300 К; ц=86,4 ГПа, у=0,29, 
а=2,867 10° М, &0=10" с", (,=0,9 эВ, х=330 МПа, с.=18 МПа, о, =2 ГПа, т=1 [20]. 

На основе результатов расчетов эволюции пластической деформации далее была вычис- 


лена временная эволюция КИН. В расчете предполагалось, что для КИН трещины имеет место 
представление [21]: 


к)=кико+Ве[ к], К®=к@+ш[к”о]. (13) 
В формуле (13) к!(®, Ки(® - см. формулу (1), а К”(Е) - поправка, учитывающая влия- 
ние пластической деформации на КИН: 
2 
к’ => [ КР(2)Ар, (ва,  и=хч. (14) 
р, 
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В формуле (14) К; - КИН трещины, соответствующий действию одной дислокации [21]. 
В некоторых расчетах выражения (13) дополняли поправкой за счет влияния области сцепления 


81 
для моды Т к КО = КИО + Ве[ К”(0 ] + К: ‚где К! о [15]. 


Результаты расчетов и их обсуждение. Проанализируем вначале результаты для модели с 
закругленной вершиной трещины. Ранее было получено [22], что распределения пластической 
деформации с учетом Т-напряжений и без них заметно различаются для всех видов ориентации 
плоскости скола, направления линии фронта трещины и плоскости легкого скольжения в кри- 
сталле. Влияние Т-напряжений при всех вариантах расчета приводит к укорочению времени дос- 
тижения равновесного распределения пластической деформации у вершины трещины в кристал- 


ле. В работе [22] для расчета с«(г,Г) использованы компоненты тензора напряжений [см. форму- 
лу (1)] при р->0 в смещенных узлах "’-+и,. В данной работе в расчетах учитывался конечный 
радиус вершины трещины (р>0) [см. формулу (1)]. Результаты расчетов показали, что форма 


пластической зоны и распределение пластической деформации в ней заметно зависят от радиуса 
вершины р, но основные особенности распределения пластической деформации (с учетом 
Т-напряжений и без них), полученные в [22], сохранились. На рис.2 показаны равновесные рас- 
пределения пластической деформации у вершины трещины моды П. Для сравнения на рис.3 при- 
ведено аналогичное распределение при р-—>0. Результаты расчетов также подтвердили выводы 
[22] о том, что наиболее наглядное представление о влиянии Т-напряжений дают временные за- 
висимости КИН и угла направления роста трещины. 
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-24 0х мкм 2,4 -2.4 о хх, мкм 2.4 


а) 6) 
Рис.2. Равновесное распределение пластической деформации без учета Т-напряжений и радиусом затупления вершины 
р = 0,12 мкм по сопряженным плоскостям скольжения. Изолинии соответствуют величине сдвиговой пластической де- 


формации (%): а- г = 0 — граница пластической области; 1 - 0,036; 2 - -0,017; 3 - 0,09; 4 - -0,07; 6- уй = 0 - грани- 


ца пластической области; 1 - -0,035; 2 - -0,07; 3 - -0,1; 4 - -0,17. Плоскости скольжения дислокаций {112}, плоскость 
трещины {100}, направление линии фронта трещины <001>, мода П 
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Рис.3. Равновесное распределение пластической деформации без учета Т-напряжений и радиусом затупления вершины 
р =0 по сопряженным плоскостям скольжения. Изолинии соответствуют величине сдвиговой пластической деформации 


х, мкм 2,4 


(%): а- Е =0 — граница пластической области; 1 - 0,18; 6- =? =0 — граница пластической области; 1 -— -—0,05; 2 - -0,1; 
3 - -0,2. Кристаллографические параметры те же, что и на рис.2 
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Исследования показали, что эти величины чувствительны также к величине радиуса вер- 
шины. Кривые 1 и 2 на рис.4 демонстрируют эффект экранирования вершины трещины для 


МОДЫ 1. 
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Рис.4. Временная зависимость коэффициентов интенсивности напряжения моды Т без учета Т-напряжений и радиусом 
затупления вершины: а- р = 0,12 мкм; 6- р=0. Кривые 1 и 2 определяют эволюцию кК.(Г) соответственно без учета 


пластической деформации и с ее учетом (экранирование вершины дислокациями). Кристаллографические параметры 
те же, что и на рис.2 
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Рис.5. Временная зависимость коэффициентов интенсивности напряжения моды П без учета Т-напряжений и радиусом 
затупления вершины: а- р = 0,12 мкм; 6- р=0. Кривые 1 и 2 определяют эволюцию к;(!) соответственно без учета 


пластической деформации и с ее учетом (антиэкранирование вершины дислокациями). Кристаллографические параметры 
те же, что и на рис.2 


Кривые 1 и 2 на рис. 5 подтверждают существование «аномальной» эволюции Ки, обна- 


руженное в [22]: при трех ориентациях плоскостей скола и линий фронта трещины при скольже- 
нии дислокаций в плоскостях легкого скольжения, образующих угол 60° с плоскостью скола, с 
учетом Т-напряжений и без них, эволюция Ки с учетом пластической релаксации дает более вы- 
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сокие значения Ки по сравнению с эволюцией Ки для хрупкой трещины. Данное явление соот- 


ветствует антиэкранированию вершины трещины. 
Кроме полученных результатов расчетов, отметим, что для данной модели (мода Г) тре- 
щины со сцеплением берегов вблизи вершины качественно повторяют особенности эволюции 
КИН за счет влияния дислокаций из пластической зоны на КИН, полученные выше. Расчеты в [15] 
выполнены в предположении «замороженности» дислокаций вблизи вершины, а в нашем случае 
происходят значительные перераспределения дислокаций у вершины трещины. 
Заключение. Затупление трещины и Т-напряжение заметно влияют на равновесное распреде- 
ление пластической деформации у вершины трещины и на временные зависимости коэффициен- 
тов интенсивности напряжения, а также приводят к уменьшению времени достижения равнове- 
сия. Подтверждено существование явления антиэкранирования дислокациями вершины затуплен- 
ной трещины, которая соответствует большим величинам Ки при учете пластической релакса- 
ции по сравнению с Ки для хрупкой трещины. Эта «аномалия» определяется плоскостью сколь- 


жения дислокаций и не исчезает при учете Т-напряжения. Результаты расчетов, учитывающих 
наличие зоны сцепления, не подтверждают оценки [15], которая получена для «замороженных» 
дислокаций. 
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РЕАЗТТС ОЕРОКМАТТОМ Е\МОГОТТОМ АТ ВЕОМТ СКАСК ТТР 1М№ СКУЗТАЕ 


О.М. КАЕРТМ$КТУ, $.\. ЗАММТКОМ 
(боиёпегп Редега! Упм№ег$Ку), 

В.М. ЗОВОЕ 

(Роп Зае Тесптса! Упмег®Ку) 


Тре р/азйс аеогтайоп ехо/иНоп аЁ те сгаск Вр т а-Ее сГу5ЕИ ипаег {пе сопаоп оЁ пе р/апе аеюгтайоп 
(тоде$ Г, П) (5 саси!ае4. Тетрога! @ётБийоп$ юг {те рас аегтайоп, ебесвуе зВеаг 5725$, 5{те55 пеп5Йу 
Гасог аге аепуеа. Тре са/сш/авоп геи! юг {пе ЫМипЕ апа аси сгасК Врз аге сотрагец. 

Кеуигогаб: сгаск Пр, р/азйс деюгтавоп, сгуяа| @!!осаНоп5, сеауаде Гасез, еазу 5[р 5узетз, °ез5 Ееп5йу 
Гасюг. 
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